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We apply the quasiparicle self-consistent GW (QSGW) method to a cubic virtual-crystal alloy 
Lai-^Ba^MnOs as a theoretical representative for colossal magnetoresistive perovskite manganites. 
The QSGW predicts it as a fully-polarized half-metallic ferromagnet for a wide range of x and lattice 
constant. Calculated density of states and dielectric functions are consistent with experiments. 
In contrast, the energies of calculated spin wave are very low in comparison with experiments. 
This is affected neither by rhombohedral deformation nor the intrinsic deficiency in the QSGW 
method. Thus we ends up with a conjecture that phonons related to the Jahn- Teller distortion 
should hybridize with spin waves more strongly than people thought until now. 

PACS numbers: 75.47.Gk 71.15.Qe 71.45.Gm 



I. INTRODUCTION 

The mixed-valent ferromagnetic perovskite 
Lai-^A^MnOs, where A is an alkaline-earth such 
as Ca, Sr or Ba, shows the colossal magnctrcsistence, 
e.g, see reviews by Tokura and Nagaosa and by 
Imada, Fujimori and Tokura 0]. As they explain, the 
colossal magnetresitance is related to the complexed 
interplay of spin, orbital and lattice degree of freedoms. 
This is interesting not only from a view of physics, 
but also of its potential applicabilities. This interplay 
can be also related to the fundamentals of high-T c 
superconductors and the multiferroic materials, which 
are now intensively being investigated [E 0] • 

In order to understand the interplay, kinds of theoret- 
ical works have been performed until now. They can be 
classified into two approaches; one is model approaches, 
and the other is the first-principle ones which are mainly 
based on the density functional theory in the local density 
approximation (LDA) or in the generalized gradient ap- 
proximation (GGA) [1 IE 0,0, El- The first-principle ap- 
proaches have an advantage that it can give energy bands 
(as quasiparticles) without any knowledge of experimen- 
tal input. Then kinds of properties are calculated based 
on the quasiparticles. However, it is well known that the 
density functional theory in the LDA (and GGA) often 
fails to predict physical properties for compounds includ- 
ing transition metals. For example, Terakura, Oguchi, 
Williams and J. Kubler [IE EH showed that the density 
functional theory in the LDA is only qualitatively correct 
for MnO and NiO. Calculated band gaps and exchange 
splittings are too small, resulting in poor a gree ment with 
optical and spin wave experiments [IE HE JUL EH ■ This 
is little improved even in the GGA. 
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As a remedy, the LDA+C/ method has been often used 
flU ]. However, it has the same shortcomings as model cal- 
culations: It can contain many parameters which are not 
determined within the theory itself, e.g, different U for 
and e g orbitals [E El and U for 0{2p) (oxygen 2p)[lj|. 
Even though there are theoretical efforts in progress to 
evaluate these U parameters in first-principle methods 
[IE |20| , we now usually have to determine these param- 
eters by hand so as to reproduce some experiments in 
practice. Then we need to check whether calculations 
with the parameters can explain other experiments or 
not. 

Many researches are performed along this line. 
Soloveyv et al. investigated LaT03(T=Ti-Cu) in the 
LDA+[/, where they tested possible ways of LDA+t/ 
in comparison with experiments. Then they concluded 
that LDA+[7 gives little differences from the results in 
the LDA in the case of LaMn03. It followed by their 
successful description of the spin-wave dispersions 0] and 
phase diagrams Q in the LDA even for i ^ 0. Ravin- 
dran et al. also showed detailed examination for LaMn03 
with full-potential calculations including spin-orbit cou- 
pling and full distortion of crystal structure Q, where 
they concluded that the density functional theory in the 
GGA worked well for LaMn0 3 . Thus, both of these 
groups reached to the same conclusion that "We can treat 
Lai_a;A x Mn03 accurately enough with the density func- 
tional theory in the LDA or in the GGA, do not need 
to use LDA+/7.". It sounds very fortunate because we 
are not bothered with difficulties about how to determine 
parameters U in the LDA+J7. However, we must check 
this conclusion carefully. For example, one of the reason 
why the GGA is accurate is based on their observation 
that their calculated imaginary part of dielectric function 
£2(w) in the GGA agrees well with an experiment Q]. 
However, this is not simply acceptable if we recall other 
cases where peaks in the calculated e2(w) are deformed 
and pulled down to lower energies when we take into ac- 
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count excitonic effects. Thus it is worth to re-examine 
the conclusion by some other methods which are better 
than those dependent on the LDA or the GGA. 

Here we re-examine the conclusion by the quasiparticle 
self-consistent GW (QSGIF) method, which is originally 
developed by Faleev, van Schilfgaarde, and Kotani [12|, 
IT3I ]. Its theoretical and methodological aspects, and how 
it works are detailed in [131 ] and references there in. They 
showed that the QSGIF method gave reasonable results 
for wide-range of materials. 

In SecHIl we explain our method. Then we give results 
and discussions in Sec lIIII In our analysis in comparison 
with experiments, calculated quasiparticle energies given 
by the QSGIF seems to be consistent with experiments. 
However, the obtained spin wave energies are about four 
times too larger than experimental values. From these 
fact, as for Lai_ x Aj;Mn03, we ends up with a conjec- 
ture that phonons related to the Jahn- Teller distortion 
should hybridize with spin waves more strongly than peo- 
ple thought until now. This is our main conclusion pre- 
sented at the end of Sec lIIH 



II. METHOD 

We first explain the QSGW method which is applied 
to calculations presented in this paper. 

The GW approximation (GWA) is a perturbation 
method. Generally speaking, we can perform GWA from 
any one-body Hamiltonian H° including non-local static 
potential V eS (r, r') as 



H 



-V 2 
2m 



+ V efi (r,r'). 



(1) 



The GWA gives the self-energy E(r, r',ui) as a functional 
of H°; the Hartree potential through the electron den- 
sity is also given as a functional of H°. Thus GWA 
defines a mapping from H° to H(oj), which is given as 

H(uj) = + V cxt + V H + E(r,r',w). Here F cxt and 
F H denote the external potential from nucleus and the 
Hartree potential symbolically. In other words, the GWA 
gives a mapping from the non-interacting Green's func- 
tion Go = l/(u> — H°) to the interacting Green' function 
G = l/(u-H(u)). 

If we have a prescription to determine H° from H(u>), 
we can close a self-consistency cycle; that is, H° —> 
H(u) -> H° -» H(u) ... (or G -> G -> G -> 
equivalently) can be repeated until converged. One of 
the simplest example of the prescription is to use H(u>) 
at the Fermi energy Ep, that is, H° = H(Ep) for 
H(uj) — > H°. In practice, we take a better choice in the 
QS GW method so as to remove the energy-dependence; 
we replace £(r, r', u>) with the static version of self-energy 
V xc (r, r'), which is written as 



(2) 



where {e^} and are eigenvalues and eigenfunction of 
H, and E(£j) = (\E , i |£(£ i )| , I'j). Re[X] means taking only 
the Hermitian part of the matrix X. With this V* c , we 
can generate a new H°, that is, it gives a procedure H° — > 
H — > H°. Thus we now have a self-consistency cycle. 
By construction, the eigenvalues of H° is in agreement 
with the pole positions of H(lo). Thus the eigenvalue is 
directly interpreted as the quasiparticle energies. This 
QS GIF method is implemented as an extension of an all- 
electron full-potential version of the GW method [2l[ as 
detailed in [13| . 

Until now they have shown that QS GIF works well for 
kinds of materials (see fl3l [22l | and references there in) . 
In [2l|, Kotani and van Schilfgaarde have shown that 
the ordinary one-shot GIF based on the LDA system- 
atically gave too small band gaps even for semiconduc- 
tors; this is confirmed by other theorists [23|, |24j. Thus 
the self-consistency is essentially required to correct such 
too small band gaps[13l |25||. Furthermore, the adequacy 
of one-shot GIF is analyzed from kinds of view points 
in (26|; e.g., it shows that the usual one-shot GIF can 
not open the band gap for Ge as shown in its Fig. 6 
(band entanglement problem). The self-consistency is 
especially important for such as transition metal oxides 
like Lai-ajAj-MnOs when reliability of the LDA and the 
GGA is questionable. We have shown that QS GIF works 
well for wide range of materials including MnO and NiO 

[n d m, m, [23, m. We observed still remaining dis- 
crepancies between the QSGIF and experimental band 
gaps, but they are systematic and may be mainly cor- 
rected by including the electron-hole correlation in the 
screened coulomb interaction IF as shown by Shishkin, 
Marsman and Kresse (23j . 

In this paper, we focus on these two objectives: 

(i) Difference of results in the QSGIF and in the LDA. 

(ii) Are results in the QSGIF consistent with experi- 
ments? If not, what can the results mean? 

For these objectives, we mainly treat the simplest cubic 
structure of the perovskite, one formula unit per cell, 
for Lai-zA^MnOa, where we set A as Ba in a virtual 
atom approximation, that is, Lai-^Ba^ is treated as a 
virtual atom with the atomic number Z = 57 — x. We 
use 6 x 6 x 6 k points in the 1st Brillouin zone (BZ) for 
integration. We also treat a rhombohedral case for x = .3 
(two formula units per cell. Its structure is taken from 
[29j ] ; angle of Mn-O-Mn is ~ 170 degree) to examine the 
effect due to the rotation of oxygen-octahedra (this is not 
a Jahn- Teller distortion). Neither phonon contributions 
nor the spin-orbit coupling are included in all presented 
calculations. 

Because of the difficulty to apply the GIF method to 
systems with localized d electrons, even the one-shot GIF 
calculations were rarely applied to Lai_ x A x Mn03 until 
now. Within our knowledge, one is by Kino et al [30|, and 
the other is by Nohara et al 3l| . Both are only within the 
atomic sphere approximation for one-body potential. In 
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contrast, our method is in a full-potential method. Thus and in the self-consistency in the QSGW 7 . 
our method here is superior to these works in this point, 

I 

III. RESULTS AND DISCUSSIONS 



Lai-xBaxMn03 x=.3 a=3.934A 




FIG. 1: (color online) The left panels are energy band at x — 0.3 in QSGW and in LDA for the lattice constant, a=3.934A(as 
used in Q) in the black solid lines. The Fermi energy Ep is at eV. Right four panels show five typical eigenvalues for different 
x, not only for a=3.934A, but also for a=3.876A. Those are shown by (red) symbols x, +, Q, □ , and ■. Correspondence to 
those in the left panels are indicated by thin (green) lines with arrows. 



In the left half of FigQJ we compare energy bands in 
the QSGW 7 and in the LDA at x = 0.3 for lattice con- 
stant a= 3.934A. The energy bands are roughly assigned 
as 0(2p), <2g, and e g bands from the bottom. In its up- 
per panels, we show labels t2 g and e g to show the assign- 
ments. The QS GW gives a band gap in the minority spin 
(T), that is, it is a half-metal, though the LDA does not. 
This enhancement of half- metallic feature in G WA is al- 
ready reported even in the one-shot GW calculations by 
Kino et al [3(| • Its implication is emphasized in a recent 
review for half-metallic ferromagnet by Katsnelson et al 
The width of the e g j band in the QSGW shows 
little difference from that in the LDA. In the QSGW 7 , 
the t2 S i band, which is hybridized with 0(2p)j, becomes 
narrower and deeper than that in the LDA. 



Right half of Fig[T] shows five typical eigenvalues as 
function of x, not only for a=3.934A, but also for 
(1 = 6. 876A. In all cases treated here, i2gT"bottom at T 
(bottom of conduction band for f ) is above Ep (Ep is 
at eV), and 0(2p)t-top at R (top of valence band for 
T) is below Ep in QSGW. This means that it becomes 
fully-polarized half metals in the QS G W (thus the mag- 
netic moment is given as 4 — x [ip,)- In contrast, the 
LDA gives a fully polarized half metal only when x = 0.5 
for a = 3.934A (i2gt-bottom is slightly above Bp). The 
eigenvalues of £2g4.-top at R are very close to that of 
0(2p)J.-top at R in the QSG1T, especially for large x. 
Though the QS GIT eigenvalues show linear dependencies 
as a function of x, the LDA does not. This is because 
the LDA has a small occupancy for the t2gT band. 
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FIG. 2: (color online) Density of states in QSGW (black solid 
line) and LDA(red dotted line) for a=3.934A. E F is at eV. 
Left panels are for minority spin, and right panels are for 
majority spin. Four panels from the top to the bottom are 
for x = 0.0, x = 0.3, x = 0.5, and for the rhombohedral cases, 
respectively. The 4/ band in QSGW is above the plotted 
region here in QSGW. 0(p z ) denotes 0(2p) orbitals along 
Mn-O-Mn. 0(p x ,p v ) are perpendicular to 0(p z ). 



Fig. [2] shows the corresponding total and partial den- 
sity of states (DOS). 0(p z ) denote 0(2p) orbitals along 
Mn-O-Mn bonding (for cr-bonding with e g orbitals). At 
first, La(4/) level is located too low in LDA, at only 
~1.5eV above Ep for x = [5|] , though the QSGW pushes 
it up to eV. At x = 0, all peak positions in 

QSGW show some disagreement with those in the LDA. 
This is due to the large difference in the occupation for 
the i 2g T band. On the other hand, for x = 0.3, 0.5, we see 
that differences of the total DOS are mainly for a peak at 
~ Ey — 2eV in |, and a peak at ~ £ , p + 2eV in f . The for- 
mer difference is related to the the t^i level. If we push 
down the LDA t2 E [ level by ~0.8eV, the occupied bands 
will be closer to those in QSGW^. The latter difference is 
related to both of the unoccupied Mn(3e£) (t2 g T and e g t). 
The QSGWpushes up £ 2g T and e g T by ~ feV, relative to 
the LDA results. The experimental position of f2g| band 
is described well in the QSGFFthan the LDA as follows. 



Angle resolved photo emission spectroscopy (ARPES) by 
Liu et al [33[ concluded that Mn(3d) band (presumably 
due to i2gl) is ~ leV deeper than the LDA result for 
La.66Cao.33Mn03. Chikamatsu et al [34j also performed 
ARPES for Lao.6Sro.4Mn03, showing that there is a flat 
dispersion around Ep — 2eV. These experiments for t2 g l 
support the results in the QSGH^. As for the positions 
of the unoccupied Mn(3e?) bands, no inverse photoemis- 
son experiments are available to identify them though we 
give some discussion below when we show 62- 

As we see above, the main difference between the 
QSGW and the LDA is interpreted as the difference of 
the exchange splitting for the i 2 g band. Roughly speak- 
ing, center of i 2 gl and i 2 gt given in the LDA is kept in the 
QSGW. Because of the larger exchange splitting, QSGW 
shows large half-metallic band gaps. In addition, the e g t 
band is pushed up. Based on the knowledge in other 
materials together with the above experiments, we think 
that the QS G W should give better description than the 
LDA. Generally speaking, LDA can introduce two types 
of errors when we identify the Kohn-Sham eigenvalues 
with the quasiparticle energies. One is the t/-type effect 
as in LDA+U. This is onsite contribution for localized 
orbitals. The other is the underestimation of the band 
gap for extended orbitals as in semiconductors (In the 
case of diatomic molecule, non-locality in the exchange 
term can distinguish bonding and anti-bonding orbitals, 
though onsite U can not). As seen in [l3j], QSGW can 
correct these two at the same time without any parame- 
ters by hand. 

FiglSl shows e 2 in comparison with the experiment [35[ 
for La^Sri-^MnOs. The LDA seemingly gives reasonable 
agreement with the experiment. For example, the peak 
position around 4eV, which is mainly due to transitions 
within the | channel, seemingly give excellent agreement 
with the experiment (upper panel). This is consistent 
with the conclusion by Ravindran et al [9j] . On the other 
hand, e 2 in the QSGW makes peak positions located at 
higher energies than the experiment by ~ 1 eV. However, 
this kind of disagreement is what we observed in other 
materials [l3j. where we identified two causes making the 
difference; (a) A little too high unoccupied quasiparti- 
cle energies in the QSGW and (b) The excitonic effect 
which is related to the correlation motion of electrons 
and holes during the polarization (we need to solve the 
Bethe-Salpeter equation). As for (a), we have an empir- 
ical procedure to estimate the error due to (a); a simple 
empirical linear mixing procedure of 80% of y xc (Eq. @) 
with 20 % of the LDA exchange-correlation practically 
worked well as shown in [22j, |28| . We have applied this to 
the case for x = 0.3 and a = 3.934. Then the t 2g J. level (x 
in Fig. [TJ) is reduced from l.I5eV to 0.93eV. This level of 
overestimation 0.2 ~ 0.3eV is ordinary for band gaps of 
semiconductors 
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If this estimation is true, the main 
cause of the disagreement should be due to (b). We think 
this is likely because we expect large excitonic effects due 
to localized electrons. At least, the disagreement in FigJ3] 
do not mean the inconsistency of the QSGW results with 
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Photon energy (eV) 

FIG. 3: (color online) Imaginary part of the dielectric func- 
tion £2 for a=3.934A. Local-field correction is neglected but it 
should be negligible as in the case ofMnOandNiO [3. Up- 
per panel is to compare calculations in QS GW and in LDA 
with an experiment for Lao.3Sro.7Mn03 [35l | by black solid 
line. In lower panel, we show results in QSGVFfor different 
x. Because of a limit in our computational method, 62 for 
< 0.5eV are not calculated (In practice, our results are with 
the wave vector q = 2^(0,0,0.04) instead of q = 0, though 
we confirmed little changes even at q — > 0.) 



SW dispersion for La-|. x Ba x MnO a=3.934 



r x r Mr R 




FIG. 4: (color online) Spin wave dispersion along F-X, F- 
M, and T-R lines for a=3.934A. Negative energy means the 
unstable modes. We also superpose the spin-wave dispersion 
in a rhombohedral case for x — 0.3 by pink lines (only T-X 
and r-R/2 in the QSGW panels). It is almost on the cubic 
case, where R/2=(0. 25, 0.25, 0.25) in cubic structure is on the 
BZ boundary of the rhombohedral structure. 



the experiment, though we need further research on it in 
future. In addition, the effect due to the virtual crystal 
approximation is unknown. We have also calculated ti 
for rhombohedral structure, resulting in very small dif- 
ferences from that for the cubic one. The lower panel in 
FigEl shows changes as a function of x. Its tendency as 
function of x (the first peak at ~5eV is shifted to lower 
energy, and the magnitude of the second peak at ~ 9eV 
is reduced for larger x) is consistent with the experiment 

Let us study the magnetic properties. As discussed 
in 0], the exchange interaction is mainly as the sum of 
the ferromagnetic contribution from the e g bands, and 
the anti- ferromagnetic one from the t2 g bands. By the 
method in [l5j , where the spin wave calculation based on 
the QSGW reproduced experimental results very well for 
MnO and NiO, we obtain the spin wave dispersions as 
shown in Fig. [4j The method is in a random-phase ap- 
proximation to satisfy (spin wave energy) — > at the 
wave vector q — > 0. In the LDA, the ferromagnetic 
ground state is stable at x = 0, but it becomes unstable 
at x > 0.3. This is consistent with the result by Soloveyv 
et al [3] , though our LDA results are a little smaller than 
those for larger x. On the other hand, we found that 
the ferromagnetic state is stable even at large x in the 
QSGW: Roughly speaking, the spin wave energies in the 
QSGW are about four times larger than experimental 
results [36[| . We also show the spin waves for the rhom- 
bohedral case in Fig. 0] (along T-X and along r-R/2), 



but they are almost on the same line in the cubic case. 
This means that the rotation of the oxygen tetrahedra 
gives little effects for its magnetic properties. In order to 
check effects of overestimation of the exchange splitting 
in the QSG , W / [13ll23l.l25l|. we use the linear mixing of the 
20% LDA exchange correlation as we already explained 
when we discuss £2, and calculate the spin- wave disper- 
sion. Then it reduces the spin- wave dispersion by ~ 11 
%, thus our conclusion here is unchanged. Our results for 
the spin-wave dispersion in the QS GlFcan be understood 
as a result of the reduction of the anti-ferromagnetic con- 
tribution of the t2 g bands because of their large exchange 
splitting. 

As a summary, our result for the spin-wave dispersions 
in the QSGW is clearly in contradiction to the exper- 
iments 136, In contrast, we have shown that 
the quasiparticle levels and £2 are reasonable and consis- 
tent with experiments. Therefore we conjecture that it 
is necessary to include the degree of freedom of phonons 
through the magnon-phonon interaction so as to resolve 
the contradiction. In fact, [3!||4(| had already suggested 
that the magnon-phonon interaction can change the spin- 
wave dispersions largely by the strong hybridization with 
phonons of the Jahn- Teller distortion. In contrast, the 
magnon-phonon interaction was supposed to play a much 
smaller role for the spin-wave dispersion in some exper- 
imental works than our result suggests: For example, 
Dai et al [37| and Ye et al [36[ claimed only the soften- 
ing around the BZ boundaries are attributed to the hy- 
bridization: Moussa claimed that the spin-wave disper- 
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sion is little affected by the magnon-phonon interaction 

In conclusion, the QS GIT gives a very different picture 
from the LDA for the physics of Lai_ x A x Mn03. The 
main difference from LDA is as for the magnitude of the 
exchange splitting for the t2 S band. It is ~2 eV larger 
than that in the LDA. QSGW predicts a large gap in 
the minority spin (i.e., it is fully polarized). Our results 
are consistent with the ARPES and the optical measure- 
ments, but not with the spin-wave measurements. We 
think that this disagreement indicates a very strong hy- 
bridization of spin wave with the Jahn- Teller type of 



phonons. It should be necessary to evaluate its effect 
based on a reliable first-principle method. 
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